Nonlocal Transformation Optics 
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We show that the powerful framework of transformation optics may be exploited for engineering 
the nonlocal response of artificial electromagnetic materials. Relying on the form-invariant prop- 
erties of coordinate-transformed Ma:xweirs equations in the spectral domain, we derive the general 
constitutive "blueprints" of transformation media yielding prescribed nonlocal field-manipulation 
effects, and provide a physically-incisive and powerful geometrical interpretation in terms of de- 
formation of the equi- frequency contours. In order to illustrate the potentials of our approach, we 
present an example of application to a wave-splitting refraction scenario, which may be implemented 
via a simple class of artificial materials. Our results provide a systematic and versatile framework 
which may open intriguing venues in dispersion engineering of artificial materials. 

PACS numbers: 42.70.-a, 42.70.Qs, 78.20.Ci, 42.25.Bs 
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Spatial dispersion, i.e., the nonlocal character of the 
electromagnetic (EM) constitutive relationships [H, @| j is 
typically regarded as a negligible effect for most natural 
media. However, there is currently a growing interest in 
its study, in view of its critical relevance in the homoge- 
nized (effective-medium) modeling of many artificial EM 
materials of practical interest jSj (based, e.g., on small 
resonant scatterers [4 m, wires [6|-|8[, layered metallo- 
dielectric composites @, [ifll, etc.), as well as in a variety 



of related effects including artificial magnetism , wave 
splitting into multiple beams 13, 12, 13|, beam tailoring 
14 1, and ultrafast nonlinear optical response [31 . If, for 
most metamaterials, spatial dispersion is seen as a nui- 
sance, counterproductive for practical applications [l6| . 
its proper tailoring and engineering may add novel de- 
grees of freedom in the wave interaction with complex 
materials 



17|. 



The transformation optics (TO) paradigm 3, 1^ has 
rapidly established itself as a very powerful and versa- 
tile approach to the systematic design of artificial mate- 
rials with assigned field-manipulation capabilities (see, 
e.g., 0, [m for recent reviews). Standard TO basi- 
cally relies on the form-invariant properties of coordinate- 
transformed Helmholtz [18] and Maxwell's equations [l9| . 
Recently, alternative approaches have been proposed, 
based, e.g., on direct field transformations '221], triple- 
space-time transformations :23] , and complex coordinate 
mapping 2J], which generalize and extend the class of 
"transformation media" (e.g., to nonreciprocal, bian- 
isotropic, single-negative, indefinite, and moving media) 



that may be obtained. Although the approach in [23] 
seems potentially capable to account for nonlocal effects, 
attention and applications have been hitherto focused on 
local transformation media. In this paper, we propose to 
apply the TO approach to develop a systematic frame- 



work for the engineering of nonlocal artificial materials, 
paving the way to novel metamaterial devices and appli- 
cations. 

Our proposed approach is based on coordinate trans- 
formations in the spectral (wavenumber) domain, where 
nonlocal constitutive relationships are most easily dealt 
with in terms of wavenumber-dependent constitutive op- 
erators. For simplicity, we start considering a distribu- 
tion of electric and magnetic sources (J' and M', respec- 
tively) radiating an electromagnetic (EM) field (E', H') 
in a vacuum auxiliary space, identified by primed coordi- 
nates r' = {x' ,y' ,z'). In the time-harmonic [exp(— icji)] 
regime, and introducing the spatial Fourier transform 



j'(k') = j G'(r') exp (-ik' • r') dr'. 



(1) 



the relevant Maxwell's curl equations can be fully alge- 
brized in the spectral (k') domain, viz., 

ik' X E'(k') = iw/ioH'(k') - M'(k'), (2a) 
ik' X H'(k') = -itj£oE'(k') -I- J'(k'), (2b) 

with eg and /^o denoting the vacuum electrical permittiv- 
ity and magnetic permeability, respectively. Throughout 
the paper, boldface symbols identify vector quantities, 
and the tilde ~ identifies spectral-domain quantities. 

We now introduce a real-valued coordinate transfor- 
mation to a new spectral domain k. 



k' = A (k) • k = F (k) , 



(3) 



with the double underline identifying a second-rank ten- 
sor operator, and the superscript ^ denoting the trans- 
pose. Similar to the spatial-domain TO, we exploit the 
form-invariant properties of Maxwell's equations in the 



2 



mapped spectral domain k [and associated, via ([T]), spa- 
tial domain r = {x, y, z)] in order to relate the corre- 
sponding fields (E, H), sources (J, M), and constitu- 
tive relationships (in terms of relative permittivity and 
permeability tensors e and jl, respectively) to those in 
vacuum [cf. ^] 

{e, h} (k) = r^{k) ■ {e', H'} [f (k)] , (4a) 

{j,M}(k)=dct^ A(k) A(k)-{j',M'} F(k) ,(4b) 



det 



1 r~ 



A(k) A(k) • A (k), 



(4c) 



with det(-) denoting the determinant, and the superfix 
denoting the inverse transpose. 
A few general considerations are in order. First, we 
note that the relationships in (U) formally resemble those 
encountered in the standard (spatial-domain) TO ap- 
proach [l^, and trivially reduce to them in the particu- 
lar case of Zmear spectral mapping [i.e., k- independent A 
in ([3])], which is fully equivalent to the local coordinate 
mapping r' = A • r. However, for a general nonlinear 
spectral mapping in ([3]), the resulting constitutive ten- 
sors in (Hcl) are always k-dependent, i.e., the associated 
constitutive relationships are nonlocal. Similar to the 
spatial-domain TO approach, the spectral field/source 
transformations in (|4a|) and (j4b|) may be used to system- 
atically design a desired response in a fictitious curved- 
coordinate spectral domain (in terms of a given nonlocal 
transformation of a reference field/source distribution in 
vacuum). Such response may be equivalently obtained 
in an actual physical space filled up by a nonlocal trans- 
formation medium whose constitutive "blueprints" are 
explicitly given by (j4c| . Restrictions on the coordinate 
mapping in ^ may be imposed so as to enforce specific 
physical properties, such as passivity and/or reciprocity. 
For instance, it can readily be verified that the Hermi- 

~T ~ * 

tian condition A (k) = A (k) yields a lossless medium, 
whereas the center-symmetry condition A (k) = A (— k) 
yields a reciprocal medium. 

In the spatial-domain TO, the choice of the coordinate 
transformation is guided by intuitive geometrical consid- 
erations essentially based on the geodesic path of light 
rays. Likewise, our nonlocal TO approach admits a ge- 
ometrical interpretation in terms of direct manipulation 
of the dispersion characteristics via deformation of the 
equi-frequency contours (EFCs). While perhaps less in- 
tuitive, such interpretation is equally insightful and pow- 
erful, as the geometrical characteristics (e.g., asymptotes, 
symmetries, inflection points, single/multi-valuedness) of 
the EFCs fully determine the kinematical (wavevector 
and velocity) properties of the wave propagation and 
reflection/refraction Figure [T] schematically illus- 

trates this interpretation with reference to an (x, z) two- 
dimensional (2-D) scenario where the EFC pertaining 



to the vacuum space is given hy k'1 + k'1 = [Fig. 
[Ija)], i.e., a circle of radius ko = uj/cq (the vacuum 
wavenumber, with cq — l/y^eo/Io denoting the corre- 
sponding speed of light). Figures [TJb) and [He) show two 
qualitative examples of transformation-medium EFCs, 



F^{k^,k,) + F^{k,,k,) 



(5) 



which, depending on whether the mapping in (|3]) is single- 
or dou&Ze-valued, may feature a moderate deformation 
[Fig. [TJb)] or the appearance of an extra branch [Fig. 
[IJc)], respectively. Our geometrical interpretation there- 
fore establishes a straightforward connection between the 
multi- valued character of the mapping and the presence 
of additional extraordinary waves. 

Moreover, for the same 2-D scenario above, assuming 
a single- valued mapping and letting $o(^) ^^id ^di^) the 
aperture distributions of a transverse (electric or mag- 
netic) field at the input {z = 0) and output (z = d) 
planes, respectively, in an unbounded material space, the 
corresponding (1-D) spatial spectra will be related via 

$d(fc.) = exp {ik,d) $o(fc.) = f{k^)^o{k,), (6) 

where the first equality arises from straightforward plane- 
wave algebra, and we have assumed 



:log 



T(fc,) 



(7) 



Equation ^ may be interpreted as an input-output rela- 
tionship of a shift-invariant linear system, in terms of the 
modulation transfer function T{kx). Within the frame- 
work of our approach, Eq. ([T]) directly defines in explicit 
form the EFC shape that is needed in order to engineer a 
desired field-transformation effect between the input and 
output planes. Accordingly, the possibly simplest spec- 
tral mapping ([3]) from the auxiliary vacuum space that 
can yield [via ([5])] such desired shape is 



F^{k^)-- 



kl + ^ log 



f{k,) , h{k,) = k. 



(8) 



The above examples highlight the intriguing perspec- 
tives of engineering the dispersion properties of the trans- 
formation medium via the vector mapping function F in 
([3]). Clearly, the practical applicability of the approach 
relies on the possibility to synthesize anisotropic, nonlo- 
cal artificial materials which, within given frequency and 
wavenumber ranges, suitably approximate the blueprints 
Icl) . While, compared with the better established 



synthesis of anisotropic, spatially-inhomogeneous mate- 
rials required by standard spatial-domain TO, this is sig- 
nificantly more challenging from the technological view- 
point, interesting nonlocal effects may still be engineered 
relying on simple artificial materials for which nonlocal 
homogenized models are available in the literature. 
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FIG. 1. (Color online) Schematic of EFCs pertaining to the auxiliary vacuum space (a) and the deformed versions obtained via 
a single-valued (b) and double-valued (c) spectral coordinate mapping. Also shown are the wavevectors and group velocities 
pertaining to a refraction scenario (details in the text). 



As an illustrative application example, we outline a 
step- wise procedure to design a nonlocal transformation- 
medium half-space so that a transversely-magnetic (TM) 
polarized (i.e., y-directed magnetic field) plane wave with 
assigned wavevector k' = [with angle Qi from the z- 
axis and associated group velocity v^(ki) = coki/|ki|, 
cf. Fig. [lla)] impinging from vacuum is split into two 
transmitted waves with prescribed directions, i.e., group 
velocity forming an angle Qt\ and Qt2-, respectively, with 
the z-axis. As schematically illustrated in Fig. [1] the 
wavevector(s) kt pertaining to the wave(s) transmitted 
into the transformation-medium half-space may be read- 
ily determined as the image(s) of the incident wavevector 
kj in the deformed ECF(s) ([S]) subject to the tangential- 
wavevector continuity htx — ^ix and to the radiation con- 
dition Re(A;tz) > 0. For a given transmitted wavevector 
kt, the corresponding group velocity (normal to the de- 
formed EFC) is given by 



^ eo/(kO-F(kt) 
|F(kO| ' 



(9) 



with J(k) = 9k'/9k denoting the Jacobian matrix of the 
transformation in ([3|), and the ± sign dictated by the 
radiation condition Re(i;gz) > 0. We first need to deter- 
mine a double-valued spectral-domain transformation (|3]) 
capable of mapping [via ([5])] the incident wavevector k^ 
into two transmitted wavevectors k^i and kt2 character- 
ized by a conserved tangential (i.e., x-) component and 
the desired group- velocity directions, i.e.. 



ktix = kt2x = kix = ko sin( 

Vgx0^tl,2) 



= tan( 



(10a) 
(10b) 



with Vg given by (jH]). A simple analytical solution to 
this functional problem may be obtained by assuming a 



variable-separated algebraic mapping of the form 

Fx{kx)^kx^/ao+a2k^, F,{k;:)^k,^JkiTh^, (11) 

with the coefficients ag, 02,^0 and 62 to be determined. 
First, by substituting (ITT|) in ([S|) [and taking into ac- 
count (|10ap ]. the above choice allows analytical closed- 
form calculation of the transmitted wavevectors kti and 
kt2, via straightforward solution of a biquadratic equa- 
tion (see f26| for details). Next, by substituting kti and 
kt2 in (jlObp [with Q and (fTTj) ] , we obtain an analytically- 
solvable system of two algebraic equations, whose solu- 
tions constrain two coefficients (say bo and 62) in (jlip . 
thereby defining a family of (infinite) coordinate transfor- 
mations which, for the given incidence conditions, yield 
the prescribed kinematical characteristics {9ti and 812) 
of the two transmitted waves (see [2^ for details). For 
instance, assuming an incidence angle 9i = 40°, and 
two transmission angles 9ti = 70° (i.e., positive refrac- 
tion) and 9t2 = —45° (i.e., negative refraction). Fig. 
^a) shows [blue (solid) curves] the double-valued EFCs 
pertaining to one such transformation (with parameters 
given in the caption). 

Note that the two seemingly free parameters in the 
transformation (oq and 02) may be in principle exploited 
to enforce additional conditions (e.g., at a different fre- 
quency). Nevertheless, we found it useful to maintain the 
flexibility endowed by such degrees of freedom in order to 
facilitate the engineering of the transformation medium 
required. Within this framework, from (|4cp . we first need 
to define the tensor operator A associated [via ([3])] to the 
vector mapping F in (ITT|) . The possibly simplest choice 
27| is the diagonal form A = diag[Fx /kx,Ayy,Fz/kz], 



where, for the assumed TM-polarization, the component 
Ayy represents a degree of freedom which may be judi- 
ciously exploited so as to simplify the structure of the 
arising transformation medium. Paralleling the spatial- 
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FIG. 2. (Color online) Examples of EFCs (a) and constitutive parameters (b,c) pertaining to a refraction scenario featuring 
the splitting of a plane wave with incidence angle di = 40° into two transmitted waves with angles 9ti = 70° and 9t2 = —45°. 
Blue (solid) curves represent the TO-based blueprints, obtained from (|lip and (|12p with ao = 0.877, a2 — — 0.0289fc^^, 
bo = 0.0934, &2 ~ — 0.0014fc^^. Red (dashed) curves pertain to the synthesized 1-D PC (unit-cell shown in the inset) with 
Ea ~ 2.752, da — O.O668A0, Eb = —2.082, db = 0.0332Ao (see 26] for details). The vertical dotted line in (a) defines the incident 
wavenumber ki^c, from which the transmitted wavevectors (marked with crosses) are determined. 



TO approach, a desirable property, which may strongly 
facilitate the scalability towards optical frequencies, is 
an effectively non-magnetic material (i.e., ^yy = 1) [28|. 
This is readily achieved by choosing Ayy — FxFz/{kxkz), 
which yields [from (I4cp and (jlip ] a uniaxial anisotropic 
medium whose relevant permittivity components assume 
a particularly simple variable-separated rational form 



1 



bo + 62^? ' 



1 



(12) 



whose behavior is shown [blue (solid) curves] in Figs. Ef^b) 
andlSl^c) for the same parameters as above. 

Interestingly, the parameterization in (jl2p closely re- 
sembles the nonlocal homogenized constitutive relation- 
ships derived in [oj for a 1-D multi-layered photonic 
crystal (PC), thereby suggesting that our TO-based 
blueprints may be approximated, at a given frequency 
and within limited spectral-wavenumber ranges, by such 
a simple artificial material. Accordingly, we consider a 
1-D PC made of alternating layers (periodically stacked 
along the x-axis) of homogeneous, isotropic materials, 
with relative permittivity Sa and Sb, and thickness da and 
db, respectively [see the inset in Fig. [UJb)]. In [§], the pa- 
rameters of the homogenized uniaxial medium [cf. ()12p ] 
were determined by matching its dispersion law with the 
McLaurin expansion of the exact Bloch-type dispersion 
law of the PC up to the fourth order in the arguments. 



Such nonlocal homogenized model establishes a simple 
analytical connection with the PC parameters, and is 
therefore very useful in our synthesis procedure. How- 
ever, in our case we developed a modified version [sim- 
ilar to (jl2p . but higher-order in k^], in order to accu- 
rately capture the parameter variations over the dynam- 
ical ranges involved (see for details). 

As a last step, we developed a semi-analytical proce- 
dure for determining the parameters of the PC approx- 
imant, based on the matching between the above non- 
local homogenized model and the TO-based blueprints 
in (see [2^ for details). For the example consid- 
ered, the above synthesis yields a PC with = 2.752, 
da = O.O668A0, Sb = -2.082, 4 = 0.0332Ao (with 
Ao = 2n/ko denoting the vacuum wavelength). Figure 
[2] compares the corresponding exact (i.e., Bloch) EFCs 
and (nonlocal homogenized) constitutive parameters [red 
(dashed) curves] with the TO-based blueprints. A sat- 
isfactory agreement is observed, especially within neigh- 
borhoods of the prescribed transmitted wavevectors kti 
and kt2 (marked with crosses), which are directly rele- 
vant to the refraction scenario of interest. 

As an independent validation of the above synthe- 
sis procedure, we carried out a finite-difference-time- 
domain (FDTD) simulation [2^ involving a finite-size 
PC slab (see [26| for details). Figure |3] shows a field 
map illustrating the splitting of an incident wide-waisted 
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FIG. 3. (Color online) FDTD-computed magnetic-field intensity map illustrating the transmission of a coUimated 

Gaussian beam in a finite-size PC (with parameters as in Fig. [2] and details in [2g|). The dashed line indicates the vacuum-PC 
interface, whereas the cyan and green thick arrows indicate the incident-beam and (prescribed) transmitted-beam directions, 
respectively. 



Gaussian beam into two transmitted beams with direc- 
tions consistent with the prescribed ones. It is worth 
pointing out that the positively-refracted beam arises 
from local effects, and may be predicted by standard 
effective-medium modeling. Conversely, as clearly visible 
in Fig. [3J the negatively-refracted beam originates from 
the excitation and coupling of surface-plasmon-polaritons 
propagating along the interfaces between the negative- 
permittivity and dielectric layers of the PC, and there- 
fore constitutes an additional extraordinary wave which 
can only be predicted by nonlocal modeling. 

For this particular example, one may argue that a semi- 
heuristic identification of the required artificial material 
and a direct optimization of its structural parameters (so 
as to approximate the desired EFCs) may have been as 
effective as the design based on the TO theory developed 
here. However, in a more general application scenario, for 
which a more complicated dispersion relation and a larger 
number of structural parameters may be desired, a direct 
optimization approach would require iterative numerical 
full-wave solutions of source-free EM problems, and may 
therefore become computationally unaffordable. Con- 
versely, the inverse process proposed here, while seem- 
ingly more involved, does not require full-wave modeling, 
and may still be carried out in a computationally-mild 
semi-analytical fashion as a parameter matching between 
the TO-based blueprints and the nonlocal homogenized 
model, similar to the above example. 

In conclusion, we have introduced and validated a 
spectral-domain-based TO framework which admits a 
physically-incisive and powerful geometrical interpreta- 
tion in terms of EFC deformation. Our approach cru- 
cially relies on the availability of nonlocal homogenized 
models, and allows the systematic synthesis of spatially- 
dispersive transformation media capable of yielding pre- 
scribed nonlocal field-transformation effects. Given the 
current research trend in metamaterial homogenization. 



with a variety of rigorous theories that allow the descrip- 
tion of the wave interaction in metamaterials in the spec- 
tral domain sO, Slj, and a correspondingly growing "li- 
brary" of nonlocal homogenized models, we expect that 
our approach may rapidly become an exciting option for 
spatial dispersion engineering in the near future. Also of 
great interest is the exploration of nonreciprocal effects, 
which our approach can naturally handle via the use of 
non center-symmetric coordinate transformations. 
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